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Abstract 

A method has been developed for obtaining equivalent linear two-body equations 
(ELTBE) for the system of many (A^) bosons using the variational principle. The 
method has been applied to the one- dimensional A^-body problem with pair-wise 
contact interactions (McGurie-Yang A^-body problem) and to the dilute Bose- 
Einstein condensation (EEC) of atoms in anisotropic harmonic traps at zero 
temperature. For both cases, it is shown that the method gives excellent results 
for large A^. 

PACS numbers: 03.75.Fi, 03.65.Db, 05.30.Jp, 67.90.+Z 



In this paper we present an approximate method of obtaining the eigenvalue solutions 
of the system of interacting N bosons using an equivalent two-body method similar to that 
used by Feshbach and Rubinov [1] for the triton {^H) three-body (A^=3) bound state. They 
[1] used both (i) the variational principle and (ii) a reduced coordinate variable ( not the 
hyperradius) to obtain an equivalent two-body equation for the three-body bound state {^H). 
For many-body problems, use of one reduced coordinate variable (the hyperradius [2]) was 
made to obtain equivalent two-body equations by keeping only a finite sum of terms of the 
hyperspherical expansion with K — Kmm {K is the global angular momentum). This method 
has been apphed to the ground state of the A^-body system composed of distinguishable 
particles or of bosons and also to nuclear bound states [3,4]. It was shown that the method 
leads to the correctly behaved nuclear bound states in the hmit of large A {A is the nucleon 
number) [4]. Recently, it has been used to describe the Bose- Einstein condensation (BEG) 
of atoms in isotropic harmonic traps [5]. 

In this paper we apply our method to the one-dimensional A^-body problem with pair- 
wise contact interactions (the McGuire-Yang A^-body problem [6, 7]) and to the dilute BEG 
of atoms in anisotropic harmonic traps at zero temperature. We show that the method gives 
excellent results for large N for both cases. 

For the A^-body system, our method for obtaining the equivalent linear two-body equation 
(ELTBE) consists of two steps. The first step is to give the A'"-body wave function ^(ri, r^, ...) 
a particular functional form 



where ^1,^2, and ^3 are known functions. We limit C's to three variables in order to obtain 
the ELTBE. The second step is to derive an equation for ^(Ci,C2,C3) by requiring that ^ 
must satisfy a variational principle 



with a subsidiary condition / ^*^dT — 1. This leads to a linear two-body equation, from 
which both eigenvalues and eigenfunctions can be obtained. The lowest eigenvalue is an 
upper bound of the lowest eigenvalue of the original N-body problem. To test our ELTBE 
method, we apply the method to the one-dimensional A'"-body problem and to the BEG of 
atoms in anisotropic traps at zero temperature in the following. 



*(ri,r2,...) «*(Ci,C2,C3) 



(1) 




(2) 
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There are only several known cases of exactly solvable three-body and four-body prob- 
lems. For N — 3 case it was shown [8] that the Faddeev equations [9] for one-dimensional 
three-body problem with pair-wise contact interactions are exactly solvable. For the one- 
dimensional N = 4 case, analytical solutions of the four-body Faddccv-Yakubovsky were 
obtained in [10]. We note that for nuclear three body systems with short-range interactions, 
the Schrodinger equation in three-dimension is reformulated into the Faddeev equations [9] 
which can be solved numerically after making partial wave expansion [11] or without par- 
tial wave expansion [12]. In the following, we consider an exactly solvable one-dimensional 
A^-body system as a test case for our method. 

For the one-dimensional N-body problem with the Hamiltonian H — —Yl^^i<Pldx} -\- 
'2cJ2i<j ^ ^'j) (with U — m = 1), the Schrodinger equation is exactly solvable. The 
bound and scattering states for this system have been found by McGuire [6] and by Yang 
[7]. For the case c < 0, there are bound states [6] for the system of N bosons with the wave 
function ^ = exp[(c/2) Y,i<j I ~ |]. The energy of this bound state is 

Eexact = -c^N{N^ - 1)/12. (3) 

The McGuire- Yang N-body problem gives a unique possibility to check the vahdity of 
various approximations made for the Schrodinger equation describing N particles interacting 

via short range potential. 

For this case, we seek for eigenfunction of in the form of 

F{p) 



p{N-2)/2 ' 



where p = jq Y^iKji^i — Xj)"^. We shall derive an equation for F{p) by requiring that ^ must 
satisfy a variational principle (2). This requirement leads to the equation 

i-$^^-^^^-n^)im-En,, (4, 

where 

Vip) = cNiN-l) ''^'f-"'\ '- (6) 

^'^^ ^ ^ V2^T{N/2-l) p ^ ' 



We note that Eq. (4) is exactly the form of the Schrodinger two-body equation in which 
^ — {F{p)/p)Yim, and a centrifugal potential energy is given by {N — 2){N — A)/{Ap^r) 
with identification of angular momentum quantum number I — N/2 — 1. Eq. (4) with the 
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Coulomb like potential V{p), Eq. (5), can be solved analytically, and we obtain for E the 
following expression [13] 

N{N-l)T{N/2-l/2) 

{N -2)T{N/2-l) ^ ' ^ ^ 

In the case of large N, using the asymptotic formulas for F function, lim|;j|^oo p(^!|!^j = 
z°'~^{l + we obtain E — —^N^ for the leading term of Eq. (6). On the other hand 

we have for large N case from Eq. (1), Eea^act = -f2^^(l + 0(^)). Therefore, for the 
McGuire-Yang N-body problem we have demonstrated that the ELTBE method, Eqs. (4) 
and (5), is a very good approximation for the case of large N (the relative error for binding 
energy is about 4.5 %). Furthermore, our approximation, Eq. (6), agrees remarkably well 
with exact value, Eq. (3), for any (the maximum value of binding energy relative error 
occurs for A'" = 3 and is about 10 %). 

Now, let us consider N identical bosonic atoms confined in a harmonic anisotropic trap 
with the following Hamiltonian 

^^-^E^^+oE ^i^l^' + ^Vi + ^lA) + E - r,.), (7) 

^''^ i=i ^ i=\ i<j 

where we assume Vint in the dilute condensate case to be the following form [14] 

Vintiri - rj) = S{ri - r^), 

with s-wave scattering length, a. 

For eigenfunction \1/ of H, we assume tlic solution for \E' has the following form 

where i^ = E£,x?, s^ = E£,»?, z' = J:1,zI 

We now derive an equation for ^{x,y,z) by requiring that ^{x,y,z) must satisfy the 
variational principle (2). This requirement leads to the equation 

= E^, (9) 

where 

h\d^ d\ , , 2 2 2n {N - 1){N - 3) ^ 1 1^ g 

^ -^^dx^^W^d^^^^^^^"" ^""'^ 4 

(10) 
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with g = t/o(27r)-3/2(r(A^/2)/r(A^/2 - l/2)fN{N - l)/2 and = M^a/m. To the best 
of our knowledge, Eqs. (9) and (10) have not been discussed in the hterature. 

For the positive scattering length case, a > 0, we look for the solution of Eq. (9) of the 
form 

= Y.<^^3k^f\^)^f\v)^k\^l (11) 

where = x^""-^^'^ e-Kv[-{x/ caf 12], $f (y) = exp[-(|///3,)V2], ^f{z) = 

gxp[— (^/7fc)^/2], and Cijk are solutions of the following equations 

l,m,n l,m,n 

where 



with g = [^^jj n, n = 2^u;m/ {2nh)Na, uj = {u!xU!yWz)K (^x = '^a;/'^, «y = '^y/^i and 

For the case of large iV, we have Xijkjmn ~ ^ii^jm^n, Hijk,imn ~ E5ii5jmSkn- For the 
ground state energy, using ^ = ^ = ^ = 0, we obtain 

E 5 2 

= (13) 



Nhuj 4 

We note that Eq. (13) is the exact ground state solution of Eq. (9) for large N. For the 
case of large N we can obtain an essentially exact expression for the ground state energy 
by neglecting the kinetic energy term in the Ginzburg-Pitaevskii-Gross (GPG) equation [16] 
(the Thomas- Fermi approximation [15]) as 

Etf 5 15 ^.1^1 

iVfii;=7<T^''"" (14) 
Comparising Eq. (13) with Eq. (14), we can see that for the case of large N the ELTBE 
method (Eqs. (9) and Eq. (10)) is a very good approximation, with a relative error of about 
8% for the binding energy. 

When the scattering length is negative, the effective interaction between atoms is at- 
tractive. It has been claimed that the BEC in free space is impossible [17] because the 



attraction makes the system tend to an ever dense phase. For 'Li, the s-wave scattering 
length is a = (—14.5 ± 0.4) A [18]. For bosons trapped in an external potential there may 
exist a metastable BEC state with a number of atoms below the critical value Ncr [19-27] . 

For the a < case, we can see that potential energy in Eq. (10), 

V(.,y,z) = ^(u,y-^u,y + c.^^) + |1 (^-1)(^'-3) (J_ + J_ + ^) _ M, (15) 

2 ' y ^ ' 2m 4 7/2 ^y^' 

for N < Ncr has a single metastable minimum which leads to the metastable BEC state. 
We note that for the case of large Ncr, the ELTBE method leads to the same Ncr as the 
variational GPG stationary theory [26]. To show this, let us consider an anisotropic trap, 
<-0x = i-Oy = uj±, ujz = Xuj±. Local minimum conditions A > 0, where A is a matrix with 
matrix elements A^j = d'^V/dxidxj, can be written for this case as 

71^/25151 -n- \^5^5l/Z2 + 0{^) < 0, (16) 

where — {2mw±/h,Ncr)z^ , 5±_ — {2muj±/hNcr)x^ , and n — 2{muj±/2T:TiY/'^Ncr | a |. Setting 
the left-hand side of Eq. (16) to zero and neglecting 0{j^) terms, we obtain the following 
equations for Ncr 

1-251 = <^i(l+8|^A2)V2, i_xHl = <5,(5^[i + (i+8|lA^)V2],^ ^ [i^(1^8|l^2)i/2]_ 

()_L d±_ d±_ 

(17) 

Eqs. (17) are exactly the same as equations for determining Ncr obtained from the variational 
GPG approach [26]. Taking the experimental values of ^Li trap parameters [28], u}^/2t: ~ 
152 Hz, and 00^/271 = 132 Hz, we obtain Ncr= 1456. This value of Ncr is consistent with 
theoretical predictions [23-27] and is in agreement with those observed in a recent experiment 
[28]. 

We note that the ELTBE method for a general anisotropic trap can be improved using 
a generalization of hyperspherical expansion 

*(ri,...r^,)= Y: (18) 

where the hyperspherical harmonics Y^l{Q,x), Y^li^y): a-^d Y^^(Q2) are eigenfunctions of 
the angular parts of the Laplace operators J^iLi w^-, Ylf=i and respectively. 
However, we do not expect a fast convergence of the expansion, Eq. (18), because of nonuni- 
formity of the convergence of the expansion of Y,i<j ^nt(ri — r^) in x, y, and z. 
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In summary, we have presented a method for obtaining an equivalent hnear two-body 
equation for the system of N bosons. We have apphed the method to the McGuire-Yang 
N-body problem and also to the dilute Bose-Einstein condensation in anisotropic harmonic 

traps at zero temperature. For both cases we have shown that the method gives excellent 
results for large N. 
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